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As rotating neutron stars slow down, the pressure and the density in the core region increase 
due to the decreasing centrifugal forces and phase transitions may occur in the center. We extract 
the analytic behavior near the critical angular velocity f2o, where the phase transitions occur in the 
center of a neutron star, and calculate the moment of inertia, angular velocity, rate of slow down, 
braking index, etc. For a first order phase transition these quantities have a characteristic behavior, 
e.g., the braking index diverges as ~ (Qq — f2) _1//2 . Observational consequences for first, second and 
other phase transitions are discussed. 
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The physical state of matter in the interiors of neutron 
stars at densities above a few times normal nuclear mat- 
ter densities is essentially unknown. Interesting phase 
transitions in nuclear matter to quark matter jj| , mixed 
phases of quark and nuclear matter (2|j|, kaon Q] or 
pion condensates [^|^], neutron and proton superfluid- 
ity [7 1, hyperonic matter crystalline nuclear mat- 
ter E], magnetized matter, etc., have been considered. 
Recently, Glendenning et al. [pi have considered the in- 
teresting question of rapidly rotating neutron stars and 
what happens as they slow down when the decreasing 
centrifugal force leads to increasing core pressures. They 
find that a drastic softening of the equation of state, e.g., 
by a phase transition to quark matter can lead to a sud- 
den contraction of the neutron star at a critical angu- 
lar velocity and shows up in a backbending moment of 
inertia as function of frequency. Here we consider an- 
other interesting phenomenon namely how the star and 
in particular its moment of inertia behaves near the crit- 
ical angular velocity where the core pressure just exceeds 
that needed to make a phase transition. We calculate 
the moment of inertia, angular velocities, braking index, 
etc. near the critical angular velocity and discuss obser- 
vational consequences for first and second order phase 
transitions. 

The general relativistic equations for slowly rotating 
stars were described by Hartle |J. We shall also make 
the standard approximation of slowly rotating stars, i.e., 
the rotational angular velocity is il 2 <C GA1/R 3 . For 
neutron stars with mass M = 1.4M Q and radius R ~ 10 
km their period should thus be larger than a few millisec- 
onds, a fact which applies to all measured pulsars insofar. 
Hartle's equations are quite elaborate to solve as they 
consist of six coupled differential equations as compared 
to the single Tolman-Oppenheimer-Volkoff equation ]Tc| ] 
in the non-rotating case. In order however to be able 
to analytically extract the qualitative behavior near the 
critical angular velocity Qq, where a phase transition oc- 
curs in the center, we will first solve the Newtonian equa- 
tions for a simple equation of state. This will allow us to 
make general predictions on properties of rotating neu- 
trons stars when phase transitions occur in the interior 



of a star. The corrections from general relativity are typ- 
ically of order GM/R ~10% for neutron stars of mass 
M ~ IAMq. The extracted analytical properties of a 
rotating star are then checked below by actually solving 
Hartle's equations numerically for a realistic equation of 
state. 

The simple Newtonian equation of motion expresses 
the balance between the pressure gradient and the grav- 
itational and centrifugal forces 



VP = -p (W + n x ft x r) 



(1) 



Here, V(r) is the gravitational potential for the deformed 
star and p the energy (^mass) density. We assume that 
friction in the (nonsuperfluid) matter insures that the 
star is uniformly rotating. Since cold neutron stars are 
barotropes, i.e., the pressure is a function of density, 
the pressure, density and effective gravitational poten- 
tial, $ = V — x r) 2 , are all constants on the same 
isobaric surfaces for a uniformly rotating star ||. We 
denote these surfaces by the effective radius, a, and for 
slowly rotating stars it is related to the distance r from 
the center and the polar angle 9 from the rotation axis 
along Cl by 

(2) 



r(a,8) = a [1 - e(a)P 2 (cos 9)] , 



where P2(cos9) is the 2nd Legendre polynomial and e(a) 
is the deformation of the star from spherical symmetry 
Inserting Eq. Q) in Eq. (Q) one obtains for small de- 
formations M the I — Newtonian hydrostatic equation 



ld_P = _ G rnia) + 2 n2c 



p da 



(3) 



where m(a) — Ait p(a')a' 2 da' is the mass contained in- 
side the mean radius a. The factor 2/3 in the centrifugal 
force arises because it only acts in two of the three di- 
rections. The equation (7 = 2) for the deformation e(a) 
is given in, e.g., Ref. flllfl . The deformation generally 
increases with decreasing density, i.e., the star is more 
deformed in its outer layers. 

In order to discuss the qualitative behavior near criti- 
cal angular velocities we first consider a simple EoS with 
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phase transitions for which Eq. (||) can be solved analyti- 
cally namely that of two incompressible fluids with a first 
order phase transition between energy density p± and P2 
(pi < P2) coexisting at a pressure Pq. The mass func- 
tion m(a) is very simple in the Newtonian limit and the 
boundary condition m(a) = M relates the star radius R 
to the radius of the dense core, i?o? &s 



R= R 



- l)R 

Pi 



1/3 



(4) 



where R = (SM/Anpi) 1 / 3 is the star radius in the absence 
of a dense core. Solving Eq. (El) gives the pressure 



P(a) =Po + \{Rl a>i(yG Pl - V) 



(5) 



for Rq < a < R. The boundary condition at the surface 
P(R) = in Eq. (g) gives 



displays the same non-analytic square root dependence 
as R (see Eq. (§). 

Let us subsequently consider a more realistic EoS for 
dense nuclear matter at high densities such as the Bcthc- 
Johnson EoS J[^| . At high densities it can be approxi- 
mated by a polytropic relation between the pressure and 
energy density: P = Kip 254 , where Ki = 0.021^^ 1154 
and po = m n 0.15fm~ 3 is normal nuclear matter mass 
density. As we are only interested in the dense core 
we will for simplicity employ this Bethe- Johnson poly- 
trope (BJP) EoS. The central density of a non-rotating 
1.4M Q mass neutron star with the BJP is ~ 3.4po- 
Furthermore, we assume that a first order phase tran- 
sition occurs at density p\ = 3.2po to a high density 
phase of density pi = 4po with a similar polytropic EoS 
P = K2p 2 54 - From the Maxwell construction the pres- 
sure is the same at the interface, Po, which determines 
K 2 = Ki(pi/ P2) 2 ' 54 ■ We now generalize Eq. (||) by in- 
cluding effects of general relativity. From Einstein's field 
equations for the metric we obtain from the I = part 



& =1 _ 2 {3 Po 



2nG Pl 



\4w Gp\R 2 



+(g - i)ff (1 - Ro/R)) (1 - Rl/R 2 )- 1 - (6) 

The phase transition occurs right at the center when 
Rq = corresponding to the critical angular velocity 
0,0 = LOo^ZitGpi where 
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t P R 
' GM 



(7) 



Generally, for any EoS the critical angular velocity de- 
pends on Pq, M, and p\ but not on p2- 

For angular velocities just below loq very little of the 
high density phase exists and Rq "C R- Expanding (|^) 
we obtain 



Ro 



2/02/pi - 1 — Ct>o 



(8) 



For oj > too the dense phase disappears and i?o = 0. 
The corresponding moment of inertia is for Rq <C R 
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where we used that the deformation is e = (5/4)lj 2 in 
the low density phase Q. However, for the qualitative 
behavior near fl only the contraction of the star radius 
R with the appearance of the dense core Rq is important 
whereas the deformations can be ignored. The contrac- 
tion is responsible for the term in the moment of inertia 
and is proportional to Rq oc (luq — lu) 3 ^ 2 near the critical 
angular velocity. Consequently, the derivative dl/dto 2 



1 dP 
~+P~da~ 



= -G 



iira-^P 



a?{l-2Gm/a) 



tt 2 t 



(10) 



where m(a) = 4ir J Q a p(a')a' 2 da' . In the centrifugal force 
term we have ignored frame dragging and other correc- 
tions of order fl 2 GM/R ~ O.lfl 2 for simplicity and since 
they have only minor effects in our case. By expanding 
the pressure, mass function and gravitational potential 
in the difference between the rotating and non-rotating 
case, Eq. ([To]) reduces to the I — part of Hartle's equa- 
tions (cf. Eq. (100) in §.) Note also that Hartle's full 
equations cannot be used in our case because the first or- 
der phase transition causes discontinuities in densities so 
that changes are not small locally. This shows up, for ex- 
ample, in the divergent thermodynamic derivate dp/dP. 

The rotating version of the Tolman-Oppenheimer- 
Volkoff equation ( |To| ) is now solved for a rotating neu- 
tron star of mass M = 1.4M Q with the BJP EoS in- 
cluding a first order phase transition. In Fig. [l] we show 
the central density, moment of inertia, braking index, 
star radius and radius of the interface (Rq) as function 
of the scaled angular velocity. It is important to note 
that Ro oc y/^lg — ri 2 for angular velocities just below 
the critical value f2o- The qualitative behavior of the 
neutron star with the BJP EoS and a first order phase 
transition is the same as for our simple analytic example 
of two incompressible fluids examined above. Generally, 
it is the finite density difference between the phases that 
is important and leads to a term in the moment of inertia 
proportional to (Oq — f2 2 ) 3 / 2 as in Eq. (||). 

The moment of inertia increases with angular velocity. 
Generally, for a first order phase transition we find for 
ft<ft (see also Eq. (g) and Fig. (0)) 



I = Jo 1 
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For the two incompressible fluids with momentum of iner- 
tia given by Eq. (|9|) , the small expansion parameters are 
ci =loI andc 2 = (5/2)uj$(p 2 /p 1 -l)/(2p 2 /p 1 - l-wg) 3 / 2 ; 
for £1 > fio the c 2 term is absent. For the BJP we find 
from Fig. (1) that c 2 ~ 0.07 ~ 2.2w§. Generally, we find 
that the coefficient c 2 is proportional to the density differ- 
ence between the two coexisting phases and to the critical 
angular velocity to the third power, c 2 ~ (P2/P1 — l) w o- 
The scaled critical angular velocity ojq can at most reach 
unity for submillisecond pulsars. 

To make contact with observation we consider the tem- 
poral behavior of angular velocities of pulsars. The pul- 
sars slow down at a rate given by the loss of rotational 
energy which we shall assume is proportional to the ro- 
tational angular velocity to some power (for dipole radi- 
ation n — 3) 



- (i/n» 

dt \ 2 



= -cn 



n+l 



(12) 



With the moment of inertia given by Eq. ( |1 1| ) the angular 
velocity will then decrease with time as 
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(n - l)t 
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yji-ij) 2 ^ + ■■■■), (13) 



for t > to- Here, the time after formation of the pul- 
sar is, using Eq. (12), related to the angular velocity as 
t ~ ioC^o/^)™" 1 and t = 7 /((n - l)Cfi£ _1 ), is the 
critical time where a phase transition occurs in the cen- 
ter. For earlier times t < to there is no dense core and 
Eq. ( |I3| ) applies when setting c 2 — The critical angu- 
lar velocity is f^o = ^oV^Pi — 6kHz for the BJP EoS, 
i.e., comparable to a millisecond binary pulsar. Applying 
these numbers to, for example, the Crab pulsar we find 
that it would have been spinning with critical angular ve- 
locity approximately a decade after the Crab supernova 
explosion. 

The braking index depends on the second derivative 
I" = dl/cPil of the moment of inertia and thus diverges 
(see Fig. (1)) as ft approaches {Iq from below 
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(14) 



For fl > flo the term with c 2 is absent. 

We now discuss possible phase transitions in interiors 
of neutron stars. The quark and nuclear matter mixed 
phase described in || has continuous pressures and den- 
sities. There are no first order phase transitions but at 
most two second order phase transitions. Namely, at a 
lower density, where quark matter first appears in nuclear 
matter, and at a very high density, where all nucleons are 
finally dissolved into quark matter. The latter does, how- 
ever, not occur if the star is gravitationally unstable at 



such high central densities. In second-order phase tran- 
sitions the pressure is a continuous function of density 
resulting in a more gentle transition and we find a con- 
tinuous braking index. This mixed phase does, however, 
not include local surface and Coulomb energies of the 
quark and nuclear matter structures. As shown in [j3||l5ll 
there can be an appreciable surface and Coulomb energy 
associated with forming these structures and if the in- 
terface tension between quark and nuclear matter is too 
large, the mixed phase is not favored energetically. The 
neutron star will then have a core of pure quark matter 
with a mantle of quark matter surrounding it and the two 
phases are coexisting by a first order phase transition. 
For a small or moderate interface tension the quarks are 
confined in droplet, rod- and plate-like structures [jsTJlsll 
as found in the inner crust of neutron stars [Q. It is 
found that normal nuclear matter exists only at moderate 
densities, p ~ 1 — 2pQ. With increasing density droplets 
of quark matter form in nuclear matter, then they merge 
into rod- and later plate-like structures. At even higher 
densities the structures invert forming plates, rods and 
droplets of nuclear matter in quark matter. Finally pure 
quark matter is formed at very high densities unless the 
star already has exceeded its maximum mass. Due to the 
finite Coulomb and surface energies associated with form- 
ing these structures, the transitions change from second 
order to first order at each topological change in struc- 
ture. If the surface tension is very small so are the surface 
and Coulomb energies of the structures and the transi- 
tions will be only weakly first order. If a Kaon condensate 
appears it may also have such structures pj|. Pion con- 
densates ||, crystalline nuclear matter ||, hyperonic or 
magnetized matter, etc. may provide other first order 
phase transitions. 

There may also be other transitions in neutron stars. 
The glitches observed in the Crab, Vela, and a few other 
pulsars are probably due to quakes occurring in solid 
structures such as the crust, superfluid vortices or pos- 
sibly the quark matter lattice in the core jL3). These 
glitches are very small Af2/f2 ~ 10~ 8 and have a char- 
acteristic healing time. In || a drastic softening of the 
equation of state by a phase transition to quark matter 
leads to a sudden contraction of the neutron star at a crit- 
ical angular velocity and shows up in a backbending mo- 
ment of inertia as function of frequency. As a result, the 
star will become unstable as it slows down, will suddenly 
decrease its moment of inertia and create a large glitch. 
Similarly, if the matter undergoes a phase transition at a 
critical temperature but is supercooled as the star cools 
down by neutrino emission, a large glitch will occur when 
the matter transforms to its normal state. If the cooling 
is continuous the temperature will decrease with star ra- 
dius and time and the phase transition boundary will 
move inwards. The two phases could, e.g., be quark- 
gluon/nuclear matter or a melted/solid phase. In the 
latter case the size of the hot (melted) matter in the 
core is slowly reduced as the temperature drops freezing 
the fluid. Melting temperatures have been estimated in 
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Jl4] , fil| for the crust and in |1 for the quark matter mixed 
phase. Depending on whether the matter contracts as it 
freezes as most terrestial metals or expands as ice, the 
cooling will separate the matter in a liquid core of lower 
or higher density respectively and a solid mantle around. 
When the very core freezes we have a similar situation 
as when the star slows down to the critical angular ve- 
locity, i.e., a first order phase transition occurs right at 
the center. Consequently, similar behavior of moment of 
inertia, angular velocities, braking index may occur as in 
Eqs. replacing Q(t) with T{t). 

In summary, we have shown that if a first order phase 
transitions is present at central densities of neutron stars, 
it will show up in moment of inertia and consequently 
also in angular velocities in a characteristic way. For 
example, the slow down of the angular velocity has a 
characteristic behavior Cl ~ C2v/l — t/to and the brak- 
ing index diverges as n(f2) ~ t^/yl — ft 2 /fig (see Eqs. 
(P^P^) ). The magnitude of the signal generally depends 
on the density difference between the two phases and 
the critical angular velocity ujq = Qq/^/2ttGpi such that 
C2 ~ (ft2/ 'pi — iVo- The observational consequences de- 
pend very much on the critical angular velocity fio, which 
depends on the equation of state employed, at which den- 
sity the phase transition occurs and the mass of the neu- 
tron star. By studying a range of angular velocities for a 
sample of different star masses the chance for encounter- 
ing a critical angular velocity increases. Eventually, one 
may be able to cover the full range of central densities 
and find all first order phase transitions up to a certain 
size determined by the experimental resolution. Since 
the size of the signal scales with Qq the transition may be 
best observed in rapidly rotating pulsars such as binary 
pulsars or pulsars recently formed in supernova explosion 
and which are rapidly slowing down. Carefully monitor- 
ing such pulsars may reveal the characteristic behavior of 
the angular velocity or braking index as described above 
which is a signal of a first order phase transition in dense 
matter. 



ACKNOWLEDGMENTS 

We would like to thank Gordon Baym, Larry McLer- 
ran, Ben Mottelson and Chris Pethick for valuable com- 
ments. 



[1] J. C. Collins and M. J. Perry, Phys. Rev. Lett. 34 (1975) 
1353; G. Baym and S. A. Chin, Phys. Lett. B62 (1976) 
241; G. F. Chapline and M. Nauenberg, Nature 264 
(1976) and Phys. Rev. D16 (1977) 450; M. B. Kislinger 
and P. D. Morley, Astrophys. J. 219 (1978) 1017; W. B. 
Fechner and P. C. Joss, Nature 274 (1978) 347; B. Freed- 
man and L. McLerran, Phys. Rev. D17 (1978) 1109; V. 
Baluni, Phys. Rev. D17 (1978) 2092; A. Rosenhauer, E. 
F. Staubo, L. P. Csernai, T. 0vergard and E. Cstgaard, 
Nucl. Phys. A540 (1992) 630 

[2] N. Glendenning, Phys. Rev. D46, 1274 (1992). 

[3] H. Heiselberg, C. J. Pethick, and E. F. Staubo, Phys. 
Rev. Lett. 70, 1355 (1993); Nucl. Phys. A566, 577c 
(1994). 

[4] D. B. Kaplan and A. E. Nelson, Phys. Lett. B291 (1986) 
57. 

[5] G. Baym and C. J. Pethick, Ann. Rev. Nucl. Sci. 25, 27 
(1975); Ann. Rev. Astron. Astrophys. 17, 415 (1979). 

[6] A. Akmal and V. R. Pandharipande, Phys. Rev. C56, 
2261 (1997). 

[7] 0. Elgar0y, L. Engvik, M. Hjorth- Jensen, and E. Osnes, 

Phys. Rev. Lett. 76, 1428 (1996). 
[8] N. K. Glendenning, S. Pei and F. Weber, Phys. Rev. Lett. 

79, 1603 (1997). 
[9] J. B. Hartle, Astrophys. J. 150, 1005 (1967). 
[10] J. R. C. Tolman, Phys. Rev. 55, 364 (1939). J. R. Op- 

penheimer and G. M. Volkoff, ibid. p. 374. 
[11] J. Tassoul, Theory of Rotating Stars, (Princeton Univ. 
Press, 1978). 

[12] H. A. Bethe and M. B. Johnson, Nucl. Phys. A230, 1 
(1974). 

[13] H. Heiselberg, in Proceedings of Strangeness and Quark 
Matter (eds: G. Vassiliadis, A. D. Panagiotou, S. Kumar, 
and J. Madsen), (World Scientific, Singapore, 1995), p. 
298. 

[14] D. G. Ravenhall, C. J. Pethick, J. R. Wilson, Phys. Rev. 

Lett. 50 (1983) 2066; C. P. Lorenz, D. G. Ravenhall, C. 

J. Pethick, Phys. Rev. Lett. 70 (1993) 379. 
[15] J. Schaffher, in proc. of "Nuclear Astrophysics", 

Hirschegg, jan. 11-17, 1998. 
[16] W. L. Slattery, G. D. Doolen and H. E. Dewitt, Phys. 

Rev. A21 (1980) 2087. 



4 



5.0 



4.0 



3.0 



2.0 



1.0 



0.0 



1 1 1 1 1 1 1 1 1 1— 

^ ^ fc 

— 

/ 

\ / 


— i — i — i — i — i — 

Pc/pO 

— R/10km 

— R (km) 
--- l/(2MR 75) 

(dl/dw 2 )/l " 

— n(co) 








::::::::: :::>A 







0.00 



0.05 



0.10 



0.15 



CO 



FIG. 1. Central density (in units of po), radii of the neutron star R and its dense core Ro, moment of inertia, its derivative 
I' /I = dl/dw 2 /I and the braking index are shown as function of the scaled angular velocity lo 2 — Q, 2 /(2wGpi). The rotating 
neutron star has mass IAMq and a Bethe- Johnson like polytropic equation of state with a first order phase transition taking 
place at density pi = 3.2po to pi = 4po- 
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